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The final stage of latex film formation was simulated by introducing donors and acceptors into the 
adjacent compartments of a cube. Homogenous and/or heterogeneous donor-acceptor distributions 
were chosen for difi^erent types of simulations. The interdiffusion of the donors and the acceptors 
within these cubes was generated using the Monte-Carlo technique. The decay of the donor intensity 
I{t) by direct energy transfer (DET) was simulated for several interdiffusion steps. Gaussian noise 
was added to the I{t) curves to obtain more realistic decay profiles. I{t) decay curves were fitted to 
the phenomenological equation to calculate the fractional mixing at each interdiffusion step. The 
reliability of the Fickian diffusion model in the case of heterogenous and homogeneous donor-acceptor 
distributions are discussed for latex film formation. 



I. INTRODUCTION 



Polymer latex particles have been utilized in a wide 
variety of applications in the coating and adhesive tech- 
nologies, biomedical field, information industry and mi- 
croelectronics. In many of these applications, e.g., coat- 
ings and adhesives, latexes form thin polymer films on a 
substrate surface. Properties (mechanical, optical, trans- 
port, etc.) of the final film should be tailor-made accord- 
ing to the application. 

Film formation from latex particles is a complicated, 
multistage phenomenon and depends strongly on the 
characteristics of colloidal particles. In general, aque- 
ous or non-aqueous dispersions of colloidal particles, with 
glass transition temperature {Tg) above the drying tem- 
perature, are named hard latex dispersion, however aque- 
ous dispersion of colloidal particles with Tg below the 
drying temperature is called soft latex dispersion. The 
term "latex film" normally refers to a film formed from 
soft particles where the forces accompanying the evapo- 
ration of water are sufficient to compress and deform the 
particles into a transparent, void-free film^'^. However, 
hard latex particles remain essentially discrete and unde- 
formed during drying process. Film formation from these 
dispersion can occur in several stages. In both cases, the 
first stage corresponds to the wet initial state. Evapora- 
tion of solvent leads to second stage in which the particles 
form a closed pack array, here if the particles are soft they 
are deformed to polyhedrans (see Figure 1). Hard latex 
however stay undeformed at this stage. Annealing of soft 
particles cause diffusion across particle-particle bound- 
aries which leads the film to a homogeneous continuous 
material. In the annealing of hard latex system, however 
deformation of particles first leads to void closure^'"* and 
then after the voids disappear, diffusion across particle- 
particle boundaries starts, i.e. the mechanical properties 
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FIG. 1: A pictorial representation of the stages of latex film 
formation from soft polymer particles, a) The latex disper- 
sion, b) The solvent evaporates leaving the particles in close 
contact, c) Deformation and packing of the particles, d) Fur- 
ther coalescence produces a mechanically rigid film. 



of hard latex films can be evolved by annealing; after all 
solvent has evaporated and all voids have disappeared. 

Transmission electron microscopy (TEM) has been 
used to examine the morphology of dried latex films^'^. 
These studies have shown that in some instances the par- 
ticle boundaries disappeared over time, but in other cases 
the boundaries persisted for months. It was suggested 
that in the former case particle boundaries were healed 
by polymer diffusion across the junction. In the last few 
years, it has become possible to study latex film forma- 
tion at the molecular level. Small-angle neutron scat- 
tering (SANS) was used to examine deuterated particles 
in a protonated matrix. It was observed that the ra- 
dius of the deuterated particle increased in time as the 
film was annealed^ and as the polymer molecules diffused 
out of the space to which they were originally confined. 
The process of interparticle polymer diffusion has been 
studied by the direct energy transfer (DET) method, 
using transient fluorescence (TRF) measurements^'^ in 
conjunction with latex particles labelled with donor and 
acceptor chromophores. The steady state fluorescence 
(SSF) method combined with DET was also used for 



2 



studying film formation from hard latex particles^"" 
An extensive review of the subject is given in reference 
14. In DET measurements distribution of donors and ac- 
ceptors are thought to be crucial i.e. it is believed that 
donors and acceptors have to be distributed randomly in 
the latex particles for the reliable TRF measurements, to 
determine the diffusion coefficients, D of polymer chains. 

Recently we have performed various experiments with 
photon transmission method using an U.V. Visible 
(UVV) spectrophotometer to study latex film forma- 
tion from PMMA and PS latexes in where void closure 
and interdiffusion processes at the junction surfaces are 
studied^^"^^. All these studies indicate that annealing 
leads to polymer diffusion and mixing as the particle 
junction heals during latex film formation. Recently, 
Monte Carlo simulation of interdiffusion and its moni- 
toring by DET during latex film formation has also been 
studied in our laboratories^^'^°. 

In this work, Monte Carlo method was used to simulate 
the final stage of film formation by introducing donors 
and acceptors into the adjacent compartments of a cube. 
Four different combinations of donor-acceptor distribu- 
tions were chosen for the different types of simulations. 
For example in the first case distribution of donors and 
acceptors in their adjacent compartments are taken as 
homogenous and gaussian respectively. In the second 
case distributions are switched from one compartment 
to the other. In the third case, both distributions are 
taken as gaussian and in the final case, distribution of 
donors and acceptors are both taken homogeneously to 
compare this case with the others. 

The interdiffusion of donors and acceptors between 
these adjacent compartments was randomly generated 
by Monte Carlo method. The decay of the donor in- 
tensity, I{t) by DET was simulated for several interdif- 
fusion steps and a gaussian noise was added to generate 
the realistic time resolved fiuorescence data. I(t) decays 
were fitted to the phenomenological equation to obtain 
the fractional mixing at each interdiffusion step. The re- 
liability of the Fickian model and the effect of heteroge- 
nous donor-acceptor distributions are discussed at the 
last stage of latex film formation process. 



II. DET AND FLUORESCENCE DECAY 
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FIG. 2: Several snapshots of the interdiffusion process be- 
tween adjacent compartments of a cube in which donors and 
acceptors are distributed homogenous and gaussian wise, a) 
K=0.0, b) K=0.3 and c) K=1.0 . 



Polymer diffusion obeys de Gennes scaling laws for 
times short compared to the tube renewal time ttr, but 
for long times it is like a random walk process (Fickian 
diffusion). In order to be able to determine whether the 
diffusion is Fickian, one must compare the experimental 
data with the results of simulations of DET with Fickian 
diffusion. 

TRF in conjunction with the DET method, monitors 
the extent of interdiffusion of donor (D) and the acceptor 
(A) labelled polymer molecules. The sample is made of 
a mixture of D and A labelled latex spheres. When this 
sample is annealed for a period of time and the donor 



fluorescence profiles are measured, each decay trace pro- 
vides a snapshot of the extent of interdiffusion^. A film 
sample after annealing was considered to be composed of 
three regions; unmixed D, unmixed A and the mixed D 
- A region. This model was first empirically introduced 
by the two component donor flourescence decay^^'^^. 

When donor dyes arc excited using a very narrow pulse 
of light, the excited donor returns to the ground state ei- 
ther by emitting a fluorescence photon or through the 
nonradiative mechanism. For a well behaved system, 
after exposing the donors with a short pulse of light. 
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FIG. 3: Several snapshots of the interdiffusion process be- 
tween adjacent compartments of a cube in which donors and 
acceptors are distributed gaussian and homogenous wise, a) 
K=0.0, b) K=0.3 and c) K=1.0 . 



FIG. 4: Several snapshots of the interdiffusion process be- 
tween adjacent compartments of a cube in which both donors 
and acceptors are distributed gaussian wise, a) K=0.0, b) 
K=0.3 and c) K=1.0 . 



the fluorescence intensity decays exponentially with time. 
However, if acceptors arc present in the vicinity of the 
excited donor, then there is a possibility of DET from 
the excited donor to the ground state acceptors. In the 
classical problem of DET, neglecting back transfer, the 
probability of the decay of the donor at rk due to the 
presence of an acceptor at is given by^^ 

Pk{t) = exp[-t/To - Wikt] (1) 

where Wik is the rate of energy transfer given by Forster 

as 
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Wik = -^K — ( — j 
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(2) 



Here Rq represents the critical Forster distance and k 
is a dimensionless parameter related to the geometry of 
interacting dipolc. If the system contains No donors and 
Na acceptors, then the donor fluorescence intensity decay 
can be derived from the equation (2) and given by^^ 
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= exp(-i/ro)-^ j nD{rk)drk 
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FIG. 5: Several snapshots of the iriterdiffusion process be- 
tween adjacent compartments of a cube in which both donors 
and acceptors are distributed homogenously. a) K=0.1, b) 
K=0.4 and c) K=1.0 . 



X exp(- J nA{ri)dri{l - exp {-Wikt))) (4) 

This equation can be used to generate donor decay pro- 
files by Monte- Carlo techniques. It is shown that the 
equation (4) reduces to a more simple form which can be 
compared to the experimental data'^. Their argument is 
summarized below for clarity. Changing to the coordi- 
nate rik = ri — rk leads to, 



— = exp(-Vro)— y nu{rk)drk 

X II / nA{rik+rk)drikex^{-Wikt){^) 

where Rg is an arbitrary upper limit. Placing a particular 
donor at the origin and assuming that the mixed and 
unmixed regions are created during interdiffusion of D 
and A, the equation (5) becomes 
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where 
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nD{rk)dr{k) 



(7) 



represent the fraction of donors in mixed and unmixed re- 
gions, respectively. The integral in equation (6) produces 
a Forster type of function^^'^^ 
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where C is proportional to acceptor concentration. Even- 
tually, one gets the following formula for the fluorescence 
intensity. 



Yij^ J '^A{ri)dri exp {-Wikt) (3) 



^ = Si exp(-VTo - C(-)V2) + B2 exp(-VTo) (9) 

-'(U) To 



Here no and ha represent the distribution functions of 
donors and acceptors. In the thermodynamic limit equa- 
tion (3) becomes^^ 



Here it is useful to define the mixing ratio K representing 
the order of mixing during interdiffusion of the donors 
and the acceptors as 



m 

1(0) 



= exp(-VTo)-^ j nD{rk)drk 



K = 



B1 + B2 



(10) 
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III. SIMULATION OF INTERDIFFUSION 

The interdiffusion of donors and acceptors between two 

adjacent compartments corresponds to the hist stage of 
latex film formation process. Here the geometry is simpli- 
fied using cubes instead of the polyhedrons, and donors 
and acceptors are randomly distributed in seperate adja- 
cent compartments in a cube. Figure 2a, 3a, 4a and 5a 
present the four types of combinations of donor-acceptor 
distributions. In Figure 2a donors and acceptors arc dis- 
tributed in the adjacent compartments in homogeneous 
and gaussian wise distributions, respectively. When these 
distributions are completely inversed, the situation is pre- 
sented in Figure 3a. Figures 4a and 5a present acceptors 
and donors both distributed in separate compartments in 
gaussian and in homogeneous wise distributions, respec- 
tively. 

Figures 2b, 3b, 4b and 5b present the picture after 
the Brownian motion of donors and acceptors generated 
for several interdiffusion steps for each combination of 
donor-acceptor pairs which are given in Figures 2a, 3a, 
4a and 5a, respectively. In each diffusion step, all the 
donors and acceptors move within a range of to 1 A° 
in any direction, but are reflected from the sides of the 
cube. After each diffusion step, the diffusion time incre- 
ments one unit. 25 x 10"^ diffusion steps were used for 
all sample simulations. The decay of donor intensity by 
DET is simulated for the configurations at the end of 
each 100 step of diffusion, therefore the diffusion process 
can be monitored quite clearly and accurately. Moreover, 
the average is taken over 10 different runs for each ini- 
tial distribution. Figures 2c, 3c, 4c and 5c present the 
final picture of the interdiffusion between two adjacent 
compartments in a cube. 

The donor decay profiles were generated using equation 
(4). The side of the cube, L, is taken as 500 A° and the 
Forster distance as 26 A°. The number of donors. No. 
and acceptors Na, are both chosen as 500. The Wik values 
for each donor-acceptor pair are obtained from equation 
(2). The parameter is chosen as 0.476, a value ap- 
propriate for immobile dyes^*^, and the donor lifetime 
is taken as 44ns. Equation (4) is then used to simulate 
the donor intensity I(t) decay profiles. /(O) = 2 x 10"^ 
is chosen and the decay profiles are obtained for a 250ns 
interval, divided into 250 channels of Ins each. Decay 
profiles at the several interdiffusion steps for both donors 
and acceptors are homogenously distributed in adjacent 
compartments are presented in Figure 6a. 

Here, one may also take into account the effect of the 
lamp profile when calculating the decay profiles^^'^°. To 
do so the decay profiles generated by the Monte Carlo 
simulation should be convolved with an experimental 
lamp profile, then the experimentally measured (j){t) is 
obtained by convolution of I{t) with the instrument re- 
sponse function L{t), as 
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FIG. 6: a) Decay curves at the several interdiffusion steps 
for both donors and acceptors are homogenous distributed in 
adjacent compartments. I) K=0.1, II) K=0.5, III) K=1.0 b) 
Noisy decay curves for the above picture. 



In this work, since we arc interested in the effect of donor- 
acceptor distributions on the interdiffusion, instead of us- 
ing experimental decay profiles we used generated decay 
profiles. This assumption is valid if one uses a delta, 6 
function light source (e.g. a very fast laser) as the lamp 
profile. In this case no convolution is needed and equa- 
tion (11) produces I{t). However, to obtain more real- 
istic decay profiles, gaussian noise can be added to the 
original decay profiles using Box, MuUer and Marsaglia^^ 
algorithm. In this algorithm, at first two gaussian num- 
bers { Ui , U2 ) between and 1 are created. Then Vi 
and V2 are calculated as shown below 



<f>{t) = I L{t)I{t - s)ds (11) 
Jo 



Vi = 2Ui-l 
V2 = 2C/2 - 1 



(12) 
(13) 
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FIG. 7: Comparison of the plots of the mixing ratio K versus 
diffusion time for different initial distributions as a) Donors 
and acceptors arc distributed homogenous and gaussian wise, 
b) Donors and acceptors are distributed gaussian and ho- 
mogenous wise, c) Both donors and acceptors axe distributed 
gaussian wise, d) Both donors and acceptors are distributed 
homogenously. 



Both Vi and V2 are distributed randomly in the range 
[— 1, 1]. S is calculated from these two numbers. 



TABLE I: Dtt^/o^ values are produced by fitting the data 
in Figure 7 to the equation (f9). The fits arc presented in 
Figure 8 for the various combinations of distributions. is 
the correlation coefficient for the fits. 



Donor 



Acceptor D7r^/a^(xlO 



Homogenous 
Gaussian 
Gaussian 



Gaussian 
Homogenous 
Gaussian 



Homogenous Homogenous 



1.14 ± 0.01 
1.37 ± 0.03 
1.41 ± 0.03 
0.66 ± 0.02 
0.91 ± 0.01 



0.995 
0.925 
0.973 
0.958 
0.991 



Each curve in Figure 7 is obtained from the average of a 
set of 10 runs. 

To test vi^hcthcr the simulated intcrdiffusion is Fickian 
or not, the planar sheet model is chosen^^. In this model 
the fraction of the diffusing substance that has diffused 
out of the planar sheet at time t is given by 



TT ^ ^ 



2exp( -2 ) (18) 



where D is the diffusion constant and a is the max- 
imum distance over which diffusion can occur. Since 
limfe^oo Ks = 1, eq.(18) can be written for n = in 
the form 



S = Vi^ + V2^ (14) 

US < 1, operation is unsuccessful and new Ui and U2 
numbers are created. If 5 > 1, Xi and X2 are calculated 
as shown below. 

M = q{^^)i (15) 

Xi = iViM)+p (16) 
X2 = {V2M)+p (17) 

Xi and X2 are mutually independent. They are gaus- 
sian numbers with an average p and standard deviation 
q. The noisy decay profiles for the homogenously dis- 
tributed donors and acceptors at several intcrdiffusion 
steps are shown in Figure 6b. 

In order to calculate the mixing ratio, K defined in 
equation (10) one should fit the generated decay profiles 
to equation (9). The decay profiles were fitted to equa- 
tion (9) using Levenberg-Marquart^^ algorithm. During 
fits the parameters C and tq are kept constant (C=l) and 
only the parameters Bi and B2 are varied. More than 
5000 decay profiles are fitted and the goodness of fitting is 
varied around < 1-5- The produced Bi and B2 values 
are used to obtain K values at each interdiffusion step. 
Figure 7 compares the plots of K versus diffusion time 
for the interdiffusions presented in Figures 2, 3, 4 and 5. 



ln(l-i^,) = -:^ (19) 

ln(l — K) values are plotted versus diffusion time in Fig- 
ure 8 and were fitted to equation (19). The fits obtained 
for all of the four combinations of distributions are shown 
in Figures 8a, 8b, 8c and 8d. The solid lines in the plots 
represent the fitting curve and the dots represent the dig- 
itized data. The diffusion constants and the correlation 
coefficients showing the goodness of fits are presented in 
Table I. 



IV. CONCLUSIONS 

Fits in Figure 8 and the values in Table I strongly 
suggest that people who work in TRF area have to be 
very careful to synthesize their latex particles which are 
labelled with the fluorescence dyes. In this work, it is 
observed that when the dye distribution is not homoge- 
nous different results in interdiffusion processes can be 
produced even the latex particles are in equal size. All 
data in Figure 8 present that interdiffusion saturates at 
the long time region. At the short time region the ini- 
tial donor-acceptor distribution is quite critical and ef- 
fects the interdiffusion (mixing ratio, K). When donors 
are distributed in gaussian wise, delay for the onset of 
interdiffusion is observed at early time region which is 
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FIG. 8: The plots of ln(l — K) versus diffusion time obtained for different combinations of distributions given in Figure 7. The 
solid lines present the fit of the data to equation (19). The slope of the solid lines produced diffusion constants which axe listed 
in Table I. The regions used for the linear fits are shown within arrows in Figures 8a, 8c and 8d. In Figure 8b all the data 
points axe used in the fit. 



obvious, since it takes some time for the donors to reach 
the acceptors to perform DET. In this case if the accep- 
tors arc distributed homogcnously, intcrdiffusion occurs 
with a single diffusion constant D, however if the accep- 
tors are distributed gaussian wise, two different intcrdif- 
fusion regimes can be observed at the intermediate time 
region. In other words, after a certain delay at early 
times, donors and acceptors meet quite fast to perform 
DET and then intcrdiffusion slows down and finally mix- 
ing saturates at longer times. 

When the donors are distributed homogcnously the de- 



lay at the short time region is quite small, especially if 
the acceptors are distributed gaussian wise, no delay is 

observed. In this case when the acceptors arc distributed 
either gaussian or homogenous wise, single intcrdiffusion 
regime is observed at intermediate time region where in 
both cases intcrdiffusion rate is similar and much smaller 
than when the donors are distributed gaussian wise (see 
Table I). 

In conclusion, if one assumes that the ideal distribu- 
tion for donors and acceptors in latex particles are both 
homogenous, then one has to expect that experimental 
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results for K should obey the picture in Figure 8d, even 
though the picture in Figure 8a looks much better i.e. 
interdifFusion starts with no delay and produces single 
interdiffusion constant. 
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